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Abstract 

In the present paper, the wavelet transform of Fractal Interpolation Function (FIF) 
is studied. The wavelet transform of FIF is obtained through two different methods. 
The first method uses the functional equation through which FIF is constructed. By 
this method, it is shown that the FIF belongs to Lipschitz class of order 5, (0 < 5 < 
1), under certain conditions on free parameters. The second method is via Fourier 
transform of FIF. This approach gives the A-regularity, (0 < A), of FIF under certain 
conditions on free parameters. Fourier transform of a FIF is also derived in this paper 
to facilitate the approach of wavelet transform of a FIF via Fourier transform. 
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1 Introduction 



Fractal Interpolation Function (FIF) was introduced by Barnsley [2] using the theory of 
Iterated Function System (IFS). It is a new and novel method to construct irregular functions 
from interpolation data. The construction of FIF depends on a functional operator which 
was investigated by Read [10] and Bajraktarevic [T]. It is now known as Read-Bajraktarevic 
operator and it defines a functional equation on FIF. In [6j [7J, the nowhere differentiable 
functions are studied as solution of functional equations. In this paper, the wavelet transform 
of FIF is studied through two different methods. The first method uses the functional 
equation of FIF to find bounds on the wavelet transform of FIF. By this method, we show 
that the FIF belongs to Lipschitz class of order 5 under certain conditions. The second 
method is via Fourier transform of FIF. This approach gives the A- regularity, (0 < A), of 
FIF under certain conditions on free parameters used in the construction of FIF. 

The present paper is organized as follows: In Section [21 a brief introduction on construction 
of FIF is given. For convenience, some points pertaining to notation and terminology are also 
noted in the same section. Following the approach of [I] , the Fourier transform of a general 
FIF is obtained in Section [3j The result given in jl] follows as a special case. Section H] 
is devoted to Wavelet transform of FIF. It is divided in two subsections - Subsection 14.11 
wherein the Wavelet transform of FIF is obtained using the recursive functional equation 
and Subsection 14. 2\ wherein the Wavelet transform of FIF is obtained using the Fourier 
transform of FIF. 

2 Construction of FIF 

In this section, a brief introduction on the construction of a Fractal Interpolation Function 
(FIF) is given. This section is based on (2j [3]. 
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Given an interpolation data {(xj, yi) G M. 2 : i = 0, 1, . . . , N}, where — oo < xq < x\ < . . . < 
xn < oo, the interval [xo, xn] is denoted by / and the smaller intervals [x k -i,x k ] are denoted 
by I k for & = 1,2, ... ,7V. The functions L k : I I k and F k : / x R -»■ R for Jfe = 1, 2, . . . , TV 
are defined as 

r / \ ,l ~ . x N x k-l ~ x x k 

L k (x) = a k x + b k = x H 

Fk(x,y) = ikV + qk(x) (2.1) 

respectively. In (12. ip . j k are free variables chosen such that \^f k \ < 1 and q k are continuous 
functions chosen such that the following conditions on F k are satisfied: 

\F k (x,y) - F k (x,y)\ < sd M ((x,y), (x,y)), s < 1 

F k (x Q ,y ) = yk-i and F k (x N ,y N ) = y k (2.2) 

where, cLm is a suitable metric equivalent to Euclidean metric in M 2 . Finally, the functions 
cjfci/xR— > I xM are defined as 

w fc (a; 3 ?/) = (L k (x),F k (x, y)). (2.3) 

It is shown in [2, [3] that ui k defined by (12.31) are contraction maps with respect to a metric 
equivalent to Euclidean metric. Consequently, 

{IxR-u; k ,k = l,2,...,N} (2.4) 

is a hyperbolic Iterated Function System (IFS) and therefore there exists an attractor A in 

N 

H(M. 2 ), the space of compact sets in IR 2 , such that A — |J cu k (A), where 0J k (A) = {cu k (x,y) : 

k=l 

(x,y) G A}. It is proved in [2], that attractor A of IFS ( 12. 4ft is graph of a continuous 
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function /:/—)• R which interpolates the given data : % = 0,1,..., N}, i.e. 

A = {(x,f(x)) : x E 1} and f(xi) = yi for z = 0, 1, . . . , N. The Fractal Interpolation 
Function (FIF) is thus defined as: 

Definition 2.1 (c.f. IHEJ/j The Fractal Interpolation Function (FIF) for the interpo- 
lation data {(xi, yi) : % — 0, 1, . . . , N} is defined as the continuous function f : I — > R, whose 
graph is the attractor of IFS {I x R; Uk, k = 1, 2, . . . , N}, where Uk are defined by (12 .31) . 

It is shown in [21 13] that the FIF / is a fixed point of Read-Bajraktarevic operator T defined 
by, 

T(g)(x) = F k (L- 1 (x),g(L^ 1 (x))), x 6 I k 

on the space (Q, dg), where the set Q is defined by Q = {g : g : I — > R is continuous, 

g(xo) = yo and g(xN) = yN} and the maximum metric dg is given by 

dg(g,g) = max|g(x) — g(x)\, g,g EQ. Hence, FIF / satisfies the functional equation 

f(x) = lk f{L-\x)) + q k {Ll\x)), x e I k and k = 1, 2, . . . , N. (2.5) 

For convenience, we assume that x — 0,x^ — 1, x k — x k _\ = for k = 1,...,N and 
2/o — Un = in the rest of our paper. So, 

1 k — 1 

a k = — and b k = — — . (2.6) 

Also, for notational convenience, the value of an empty product is taken to be unity and the 
value of an empty sum is taken to be zero. 



1 



3 Fourier Transform of a FIF 



In this section, the Fourier transform of a FIF is obtained for a general q^. The result given 
in [1] follows special case. 

Let / be a FIF obtained from an interpolation data. Extend / to R by defining f{x) = if 
x ^ I and hence / £ Li(R). The Fourier transform of / is given by the following theorem: 

Theorem 3.1 Let f be a FIF obtained from the interpolation data {(xi,yi) £ R 2 : % = 
0, 1, . . . , N}, where = x < X\ < . . . < x N = 1, x fc — Xk-i = -h for k = 1, 2, . . . , N and 
Ho — Un = 0. Extend f to R by defining f(x) = if x £ I = [xq, Xn]. The Fourier transform 
of f is given by 



N 





j=l k 1 ,k 2 ,...k j =l 



I 



where, /;/,../,,.../,, 



+ 



•2-1 

N 2 



+ ...+ 
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Proof The Fourier transform of / is 




Using (12. 5 p in the above equation, 
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Substituting x = Lk(x) 



x+k-l 
N ' 



we have 



k=l j 



Consider any j = 1,2,.... Replace u by jfe in (13.21) . Then 



i N 

!(-) = -y 

k=l 



-iuj(k-l)/NJ +1 



Ikf 



UJ 



i+i 



q k (x)e- iuJx/NJ+1 dx 



(3.2) 



Therefore, by induction, for any n — 1,2, 



1 N 

= E 



fei,fc2,...,fc,i = l 

AT 



n f ( 

VJV" 



7 fcl 7fc 2 • • . 7fc . .e-***!^.-^ / 9fe .(x)e-^/ iV ^x. (3.3) 



E jyj E 

j=i fei,fc 2 ,...,%=i 



With the argument similar to 0], the Fourier transform of / given by (13. ip is obtained 
from (13. 3 p as n — > oo. ■ 

Corollary 3.1 If q k ,k = 1,...,N in (12. ip are polynomials of degree i.e, qj.{%) = 
Ck, r x r ,Ck imk 7^ 0, the Fourier transform of FIF f is given by 

r=0 



oo TV 

/(<*>) = E E 7fci7fca • • • 7 ^-i e 
j=i fei,fc 2 ,...,fe 3 -=i 



3 X 



or 



x E%- 

r=l 

(-z) r+1 iV r V! 



r 7Vi zr(r - l)N 2 i r(r - l)(r - 2)iV 3 ^' 



+ . . . 



r+l 



V7V , (-z)^jVH r! 



r+l 



(3.4) 
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Proof Using qk{x) = c^ r x r for k = 1, 2, . . . , N in (13. ip and integrating, we have 

r=0 



= E 

J=l fcl,fc2,...,fe;,*=l 

E c fcj,r 
r=0 



7fci7fc 2 • • -7%-ie l£jpfe i' fe 2.-.%-i x 



x 



ziV^ A^'r iN^r{r - 1) iV 4 J'r(r - l)(r - 2) 



+ 



Co" : 



Co' 



r+1 



—iujkj 



(_A\r+l ]\f(r+l)j r \ 



Co' 



r+1 



(3.5) 



5y f j2.2Q . the constants Ck, r , k = 1,2, . . . , N satisfy the following conditions: 



ci,o = 



• C k,r = Cfc+1,0, fc = 1, . . . , N - 1 

r=0 

r=0 

Using the above conditions in f 13 . 5 1) . £/ie Fourier transform of FIF f as given by (13 ,4p is 
obtained. ■ 



Remark 3.1 m (12. ip are linear polynomials i.e, qk(x) = CkX+dk, the Fourier transform 
of the FIF given in j^j as 



oo ( N 

= ^ E - 1) ^ E ^7 fcl • • • 7*^ 

3=1 ^fci,fc 2 ,...,fcj=l 



(3-6) 



follows from above corollary. 
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4 Wavelet Transform of FIF 



In this section, the wavelet transform of a FIF is studied through two methods - first using 
the recursive functional equation given by (12. 5p and second using the Fourier transform of 
FIF given by (pTIj) . 

The wavelet transform W^g of a function g G Li(R) with respect to a suitable wavelet ip is 
a function over the half-plane H = {(s, t), s, t G R, s > 0} defined as follows: 

W^(M) = - / (~7~) dx ' 

4.1 Wavelet Transform of FIF via functional equation 

For k = 1, 2, . . . , N, assume qu G Lip 5, < 8 < 1; i.e., for some constant i\" > 0, |gfc(x) — 
9fc(y)| < A"|x - for x,y G E. Then, for k = 1,2,..., N, q k o L^ 1 : 4 -> R is a 
function defined on the compact interval If. and q k o L^ 1 G Lip 5; in fact for if* = KN 5 , 
\q k o L^ x (x) - q k o ^(y)] < K*\x - y\ 5 for x,y G R. 

Let ^ be a wavelet such that -0 G Li(R), J ip{x)dx = and defined as = x s ip(x) is 
also in Li(R). We also choose the wavelet such that the following conditions are satisfied: 



(i) ip is real and supp-0 C R + . 

(ii) For some r > 0, ip{uj) = u r + 0(a> r+1 ), a; — > + . 
(in) For each p > 0,i[)(uj) = 0{u~ p ), uj — > oo. 



LI] 



It is well known that wavelet transform of any bounded function g with respect to ip G L\ (R) 
is bounded. Since q k ° L^ 1 : I k M. are functions defined on the compact interval l k ,k = 
1,2, . . . , N, q k o L^ 1 for k — 1, 2, ... ,n are bounded and therefore each W 7 ^^ o L^" 1 ) is also 
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bounded. For a bounded function g, if g(s) = 0(\s\ 5 ) as s — > then g(s) = 0(|s| 5 ) for all s. 
Then, from [8], it is observed that the wavelet transform of q k ° Z/^T 1 satisfies 



\W4q k oL- 1 )(s,t)\ = 0(s 5 ) ii 5<r 
\W4q k oL^)(s,t)\=0(s m ) if <5>r. 

Now, for (s, t) G H, the wavelet transform of FIF is given as 



(4.2) 



/ fe=1 4 



Using ( 12 .5p in the above equation, we get, for (s,£) G if, 
W^/M) = \Y,J [l k }{L- k \x)) + q k {L- k \x))]^ {^-^ dx 

=- s it{^J \ K ) dx+ 1 (^) dx 

=t\*hl ( x ~ (m N~s k ~ l) ) <fa+ J / (^) * 

I Ik 

= V k W i>f( Ns > Nt-(k- 1)) + W4q k o L^)(s, t)\. 



k=l 



By induction, for any n = 1, 2, . . ., we have, for (s, £) G if, 



V n 

= E ^7*. • • • 7fe„ W^/(JV» aj iV n t - Nn - j (kj ' 1)) 

n AT j'—l 

+ E E TfeTfe-"^-^^^ (4.3) 

3=1 ki,k2,...,kj=l p=l 
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Defining flj = N s \ , jj\ and Q = max{f2.,- : j = 1,...,N} as in [5], an upper bound on 
\W^f(s, t)\ as s — > is obtained in the following theorem: 



Theorem 4.1 Let f be a FIF obtained from the interpolation data {(xj,?/j) G M 2 : i = 
0, 1, . . . , iV} ; where = xq < x\ < . . . < xn = 1, x^ — Xk-i = jj for k = 1, 2, . . . , N and 
Ho — Vn = 0. For this, the free parameters 7& and functions qk,k = 1,2, . . . , N in (12.11) 
satisfy, for some constant K > and < 5 < 1, < j^i+t an d \lk{x) ~ Qk(y)\ < K\x — y\ 5 
for x, y G M i.e., G Lip 5. Extend f to R by defining f(x) = if x (jL I = [xq,x^]. The 
wavelet if; is chosen such that ifi G Li(R), J ip(x)dx = 0, <f> defined as 4>{x) = (x s ip(x)) is 
also in Li(M) and the conditions (i),(ii) and (Hi) in (14. ip are satisfied for some r,p > 0. 
Then the following holds: 

(a) If5< r, for (s,t) G H, \W^f(s,t)\ < £f^| 5 , 

(b) f belongs to Lipschitz class of order 5 if 5 < r. 

Proof (a) From (@~2D, |W v ,(g fc oL fc ' 1 )(s, t)\ < K*s & if 5 < r. Substituting this bound in f l4~6D . 

we have, 

N n 

mf(s,t)\ < Yl w 1^1 ■■■ 1^1 mf(N n s,N n t-Y,N n ^(k 3 -i))\ 

ki,k2,...,k n =l j=l 
n N N 

+ E E w iT fc j ••• i7 fc ,_j J2 K *( Nj ~ ls ) s 

j=l ki,k2,---,kj—i=l kj=l 

N n 

< E iTlbil ••• ItJ |^/(iV" S ,^t-^iV^(A;,-l))| 

ki,k2,...,k n =l j=l 
n 

+ ^Ns 5 ^^)^ 1 - ( 4 - 4 ) 

Since |7fc| < a < for all k — 1, . . . , N, NQ < 1. Also, the conditions on ip tells 
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us that |W^,/| is a bounded function. Hence, as n — > oo, |W^,/(s,t)| < ^Mi \ s \ 5 f° r a ^ 
values of (s, t). 

(b) By Theorem 2.1.1 of Chapter 4 in [9], |W,/,/(s,£)| = 0(s 5 ) implies / belongs to Lipschitz 
class of order S. 



4.2 Wavelet Transform of FIF via Fourier transform 

It is well known [8] that the wavelet transform of a function g G L±(K) with respect to a 
wavelet ip G Li(R) is also obtained by the following expression: 



W^g(s,t) = — I g(u)ifj(su)e u "du. 



Let qki k — 1, 2, . . . , N in (12. ip be polynomials of degree i.e qk(x) = Ck, r x r - Choose 

r=0 

if) such that ip(oj) = M+ 7 with M — 1 > m = max{mfc, fc = 1,2, ... , iV}. Then, for 
(s, t) G H, the wavelet transform of FIF is given by 



Using (13 .4p in the above equation, we get, 



°? r oo n 

= £ 7* 

o I i =1 fci,fc 2 ,...,fcj=i 



"»k. 



or 



x £c fcj 

r=l 

(_ i )r+l jV ri r | 



7fe • ■■Ik.^e 
rW ir{r - 1)N 2 ' r(r - l)(r - 2)i\^' 



+ 



r+l 
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Hence, for (s,t) G H, 

mf(s,t)\ 



M 

< — 



OO TV T 

oo N 



E E I7fcill7fe 2 l • • • l7Ay_i|e" 

A — 1 f„ . 7„- U 1 



j=l ki,k2,...,kj=l 

E \ c ^,r\ \ rN j u M ~ 2 + r(r - l)N 2j oo M - 3 + ... + 2N rj r\ Z"^ 1 \du. 



o 



X 

r=l 



Fix < s < oo. The above integral is divided into \ W t pf(s,t)\i and \W^f(s,i)\2, where 

|MVM)li 

M l oo N 



S 

2^ 



E E l7 fel ||7 fa |---l7^Je— X 

„' — 1 U 7, 1 



j=l ki,k2,...,kj=l 

E rN j u M - 2 + r(r - l)N 2] uj M - z + r(r - l)(r - 2)N 3j co M - 4 + . 



o 



X 

r=l 



2N n r\ u M ~ r ' v }duj (4.5) 



and 



|WV/(a,*)| 2 

7 E E i^ii7 

J A — I U 7-_ 7„ 1 



S 

2^ 



OO AT 

oo N 



k 2 \ 

3=1 ki,k2,---,kj=l 



X 



17^-Je— x 

E \ c ki,A \ rN j u M ~ 2 + r(r - l)N 2j L0 M ~ 3 + r(r - l)(r - 2)N' ij u M ~ A + ... 



r=l 

.M-r-1 



+ 2N r ' J r\ u M - r - l )du (4.6) 
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Using the regularity notation X a as in [9] and from ( 14. 5 p and (14. 6p . the regularity of FIF is 
obtained in the following theorem. 

Theorem 4.2 Let f be a FIF obtained from the interpolation data {(xi,yi) G M 2 : i = 
0,1,..., N}, where = xq < x± < . . . < x^- = 1, — Xk-i = for k = 1, 2, . . . , N and 
Ho = Un = 0. Here, for k = 1, 2, . . . , N, the functions qu in (12.11) are polynomials of degree 
rrik and the free parameters jj. satisfy < N ^+i , where m = maxjmfc, k = 1,2,..., N}. 
Extend f to R by defining f(x) =0 if x G" I = [x®, xn]- Let if) be such that i/j(uj) = \co\+ e~ w 
with M — 1 > m. Then the following holds: 

(a) \W^f{s,t)\ = o(s) as s ->■ 0, 

(b) f is of regularity M — m i.e. f G A(K), A = M — m. 
Proof (a) It is observed that, for p = 2, 3, . . . , r + 1, 

s oo 

s M J u M - p e~ su) du = and s M J u M - p e~ SU] du = s p ~\M - p)\. (4.7) 

s 

Using (HID in (03]) and (Ojl . we have \W^ 3 f{s,t)\ 1 = and 

. oo AT 

|W^/(s,t)| 2 < — ^ ^ l7*illTtal ---iTfcj-Jx 
i=i fci,fe2,...,fc,'=i 



m fe 

X 

r=l 



J] |c fcj , r ||riV J (M - 2)! s M - x + r(r - l)iV 2j (M - 3)! s M " 2 + . . . 
+ 2iV rj r!(M-r- 1)! s M " r 

Now, as s -> oo, \W^f(s,t)\ = \W^f(s,t)\i = and as s -> 0, |W^/(s,£)| - 
\Wt/}f(s,t)\ 2 < Cs M ~ m for a suitable constant C. Hence, \W^f(s,t) \ = o(s) as s — > 0. 
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(b) By Theorem 2.1.1 of Chapter 4 in [9], \W^f(s,t)\ = o(s) implies / is of regularity M — m 
i.e. / e A(R), X — M — m. 

■ 
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